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1. Introduction 

Top-quark pair production is one of the key measurements performed at the LHC. The 
detailed analysis of the top quark properties will contribute to a better understanding of 
the origin of particle masses and electroweak S3niimetry breaking, and it might also give 
hints on the physics that lies beyond the Standard Model. With the large number of top 
quarks being produced at the LHC [1], the study of their properties is becoming precision 
physics. Some observables, such as the total cross section for tt production, are expected 
to be measured with an accuracy of the order of five percent. In order to match these very 
precise experimental measurements with equally accurate theoretical predictions, next-to- 
next-to-leading order (NNLO) corrections must be considered. 

Although a full NNLO calculation of the tt cross section including all required par- 
tonic channels is so far missing, a rapidly increasing number of pieces and intermediate 
results have become available recently [2-13]^. Most notably, the inclusive total hadronic 
tt production cross section induced by the all-fermion partonic processes has been com- 
puted [4,11]. 

^An up-to-date review of these intermediate results can be found in [14]. 
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At NNLO, the calculation of any observable receives three classes of contributions: 
double real da^^, mixed real- virtual, da^^, and double virtual contributions da^^ . For 
an m-jet observable, these contributions contain respectively (m + 2), (m + 1) and (m) 
partons in the final state. While the latter contribution is already in an m-jet final state 
configuration, the double real and real virtual classes of partonic channels contribute to 
the m-jet observable at NNLO if the partons present in these channels are theoretically 
unresolved (soft/coUinear) or are experimentally unresolved, i.e clustered to form an m-jet 
final state by a given jet algorithm. In addition, for hadronic observables, mass factorization 
counterterms which will enter at the (m + 1) and m-parton level have to be taken into 
account. The integration of the matrix elements with real radiated particles over the soft 
and/or collincar regions of phase space yields infrared divergencies. Therefore, in order 
to evaluate hadronic observables including higher orders, a process independent procedure 
which enables the extraction and cancellation of those infrared poles amongst the different 
partonic channels, needs to be applied. 

Subtraction methods explicitly constructing analytically integrable infrared subtrac- 
tion terms which reproduce the behaviour of the full matrix elements in their unresolved 
limits are well known solutions to this problem [15-29], [30-38], [2,39-42]. For QCD ob- 
servables involving massive fcrmions, fewer subtraction terms are needed, since the real 
radiation amplitudes are singular in fewer regions of phase space than their masslcss coun- 
terparts. Indeed for those observables, QCD radiation emitted off a massive leg can only 
lead to soft singularities. Strict collinear divergencies cannot occur, since they are regu- 
lated by the mass of the fermion. In this context, collinear limits are to be replaced by 
their massive analogues: the quasi-collinear limits [17,39]. In these limits, when integrated 
over the appropriate phase space, the real radiation matrix element is not divergent but 
develops terms of the form log{mQ/Q^), with Q being the hard scattering scale in the 
problem under consideration. Since those logarithms are not enhanced in the context of 
tt production at the LHC, they have not been taken into account in the extension of the 
antenna subtraction method developed in [2] . There, subtraction terms capturing only the 
strict infrared behaviour of the matrix element squared were derived. Even though, the 
infrared structure of the matrix elements for processes involving massive fermions is sim- 
pler, the kinematics and the integration of the subtraction terms become more involved, 
given the fact that the finite parton masses introduce a new scale into the problem under 
consideration. 

At the NLO level, two different subtraction methods have been extended to deal with 
massive final-state fermions: the dipole formalism [17, 26] and the antenna subtraction 
method [39,42]. The latter has been further extended to NNLO and employed to construct 
subtraction terms for the double real corrections to heavy quark pair production in [2,40]. 
We shall follow this second subtraction framework in this paper. 

Based on the universal factorisation properties of QCD colour-ordered amplitudes 
squared in their infrared limits, antenna subtraction [2, 30 42] constitutes a framework 
to construct subtraction terms that approximate the double real and mixed-real virtual 
corrections in their infrared limits. Within this formalism, the subtraction terms are built 
as products of antenna functions and reduced matrix elements squared with remapped 
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momenta. 

The antenna functions, which are derived from physical colour-ordered matrix ele- 
ments, capture all unresolved radiation emitted between two hard radiators and reduce to 
the well-known universal infrared factors (coUinear splitting functions or soft eikonal fac- 
tors) in the appropriate limits. Depending on whether the two hard radiators are located in 
the initial or in the final state, three types of antennae are needed: final- final (f-f), initial- 
final (i-f) and initial-initial (i-i). The initial-final and initial- initial antennae are obtained 
from their final-final counterparts by crossing one or two final state partons respectively to 
the initial state. Furthermore, by the flavour of the hard radiators we distinguish between 
quark-antiquark, quark-gluon and gluon-gluon antennae. While at NLO the subtraction of 
the unresolved limits of the real corrections only requires tree-level three-parton antennae, 
at NNLO we need tree-level four-parton antenna functions as well as products of two tree- 
level three-parton antennae for the double real contributions, and one-loop three-parton 
antennae for the mixed real-virtual contributions. 

In addition to the antenna functions and the phase space mappings needed for the 
reduced matrix elements, the method uses a Lorentz invariant factorisation of the full phase 
space into an antenna phase space, and a reduced phase space with remapped momenta. 
This factorisation is needed in a different form in each configuration (f-f, i-f, i-i), and it 
enables the analytic integration of the antenna functions over the antenna phase space. The 
integrated subtraction terms therefore contain these integrated antennae, and the reduced 
matrix elements and phase space, both defined in terms of remapped momenta, which are 
left unintegrated. This integrated subtraction term can then be combined with the virtual 
contributions and mass factorisation counterterms, achieving thus an analytic cancelation 
of the infrared poles. 

In the massless case, the phase space factorisations have been derived in [31,34]. All 
the massless final-final antennae were integrated in [34], the initial-final ones in [31,32], and 
the initial-initial antennae in [30,31,33,37]. For the antennae involving massive partons, 
the three-parton tree-level antennae, which are the only ingredients needed at NLO, have 
been integrated in [39,42]. At NNLO, so far only one antenna function with massive 
particles has been integrated [41]: the final-final four-parton antenna involving a massive 
and a massless quark-antiquark pair and derived from the physical process 7* QQqq. 

It is the purpose of this paper to evaluate the integrated forms of the new four-parton 
antenna functions employed in the double real subtraction terms for tt production due 
to the partonic processes qq — >■ ttq'q' and gg — > ttqq [2,40]. These partonic processes 
constitute the Np double real contributions to the tt hadronic cross section evaluated at 
NNLO. The three specific antennae have a massless initial-state parton and a massive 
final-state fcrmion as their hard radiators, and they are therefore initial-final antennae. 

The paper is organised as follows. In section 2, we briefly recall the infrared structure 
of double real radiation corrections to hadronic jet observables involving massive final state 
fermions. We specify in particular how the initial-final four-parton massive antennae enter 
in the construction of the subtraction terms and present the phase space factorisation that 
enables the integration of these antennae. In section 3 we recall the precise definitions 
of these specific three initial-final four-parton antennae whose integration constitute the 
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core result of this paper. Section 4 describes how the integration over the massive initial- 
final antenna phase space of these four-parton antenna functions is performed. As the 
expressions of the integrated antennae are too lengthy to be presented here, we shall give 
the corresponding pole parts in section 4, and give the full result in a Mathematica file 
attached separately to this paper. Finally, section 5 contains our conclusions. An appendix 
containing the unintegrated forms of the initial-final four-parton antennae is included too. 

2. Initial- final antenna subtraction with massive fermions at NNLO 

At the partonic level, the double real emission contributions to an m-jet cross section 
involving a pair of massive fermions (QQ) read 

m 

X-^ \Mm+4iPQ,PQ,P5, ■ ■ ■ ,Pm;Pl,P2)\'^ J^'^'^\pq,Pq,P5, . . . ,Pm) (2.1) 

= J^NNLO XI d^m+2(p3, ■ ■ ■ , Pm+i; Pi , P2) 
m 

X-^ \Mm+iijPZ, ■ ■ ■ ,Pm+4.,;Pl,P2)\^ J^^'^Hps , ■ ■ ■ , Pm+4) , (2.2) 

where the last line is obtained by relabelling all final state partons. In eq.(2.2), Sm is a 
symmetry factor for identical massless partons in the final state, |A^m+4p denotes a colour- 
ordered tree-level matrix element squared with m -|- 2 final-state partons (two of which are 
massive) and two arc initial-state partons, and denotes the sum over all the possible 
colour orderings. The next-to-next-to leading order normalisation factor J^^^^o includes 
all QCD-independent factors as well as the dependence on the renormalised QCD coupling 
constant a^. d^rn+2 denotes the phase space for an m -|- 2-parton final state containing 
m massless and two massive partons with total four-momentum Pi + P2- Finally, the jet 
function denoted by J^^"^^ ensures that out of m massless partons and a QQ pair, an 
observable with a pair of heavy quark jets in addition to (m — 2) jets, is built. 

In principle, only the leading colour pieces of the double real corrections are accounted 
for in eq.(2.2), since the subleading colour contributions involve in general interferences 
between sub-amplitudes with different colour orderings. To keep the notation simpler, 
however, we denote these interference terms also as |A^m+4p- 

The NNLO contributions given in eq. (2.2) contain infrared singularities which arise 
when one or two final state partons are unresolved (soft or collinear). The phase space in- 
tegration can only be carried out numerically after those singularities have been extracted 
using subtraction terms that approach the real radiation matrix elements in all their unre- 
solved limits. Depending on the colour connection between the unresolved partons in each 
sub-amplitude squared in eq.(2.2), the antenna subtraction method distinguishes between 
the following configurations [34, 38] 

(a) One unresolved parton but the experimental observable selects only m jets. 

(b) Two colour-connected unresolved partons (colour-connected). 
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(c) Two unresolved partons that are not colour- connected but share a common radiator 
(almost colour- unconnected) . 



(d) Two unresolved partons that are well separated from each other in the colour chain 
(colour-unconnected) . 

(e) Compensation terms for the over-subtraction of large angle soft emission. 

For each of the configurations listed above, the antenna subtraction terms have a 
distinct characteristic structure in terms of required antennae functions, which is valid 
for final-final, initial-final and initial-initial configurations. These characteristic structures, 
first derived in the massless case in [32,34,38], are unaltered by the presence of non- 
vanishing parton masses [2] , although some of the antenna functions needed in this latter 
case are massive. 

With the exceptions of subtraction terms of type (6) and (e), all the remaining configu- 
rations are constructed with products of the NLO three-parton tree-level antenna functions. 
Configuration (e) only arises in partonic processes involving at least three gluons, and its 
treatment is left to be discussed elsewhere. We shall here focus only on configuration (6) 
in the presence of massive final-state fermions. In particular, we shall concentrate on those 
cases where the two unresolved partons are colour-connected to one of the massive final- 
state particles and to one of the (massless) incoming partons. The form of these specific 
type of initial-final subtraction terms with massive antennae has been derived in [2, 40] 
in the context of the evaluation of the double real radiation corrections to top quark pair 
production at the LHC for different partonic processes. We shall now recall the form of 
these subtraction terms. 

When two unresolved partons j and k are adjacent and colour-connected to two hard 
radiators labelled i (massless) and I (massive), the initial-final subtraction term da^^^Q 
related to the double real contribution da^^j^Q given in eq.(2.2) reads [2] 

^^NNLO = ^NNLO y]d$,n+2(p3, • • • ,Pm+A\Pi,P2)-^^ 

El Y^ Y^ \ 

jk 

X \Mm{P3,---,PL,---,Pm+4;XiPi,P2)\^ J^\p3,---,PL,---,Pm+4) ■ (2.3) 



This subtraction term involves: the phase space for the production of (m -|- 2) partons, 
d$m_|_2, with two of them being massive, the colour-ordered reduced (m -|- 2)-parton ampli- 
tude squared |Alm+2p (with two partons less than the original matrix element squared), 
the jet function , the initial-final four-parton antennae Xijki and products of three- 
parton antennae in final-final and initial-final configurations. The initial-final antenna 
functions are defined by crossing one massless partons in the corresponding final-final an- 
tennae. Since the three and four-parton final-final antennae are denoted by Xj^j^ and X^^jj^-i 
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respectively, the corresponding three and four-parton initial-final antennae are denoted by 
Xl,, and ^0^.,, 

By construction, the subtraction term in eq.(2.3) contains all colour-connected double 
unresolved limits of the (m + 4)-parton colour-ordered matrix elements squared. Those 
double unresolved limits are captured by the four-parton antennae Xij^h which are one 

of the new ingredients appearing in the NNLO double real subtraction terms. The hard 
radiators are the initial-state parton i and the (possibly massive) final-state particle and 
the unresolved particles are j and k. These four-parton antennae, however, are generally 
also singular in single unresolved limits of j oi k, where they do not reproduce any physical 
singularity of the matrix element. In order to ensure that these subtraction terms are only 
active in the double unresolved regions that they arc aimed at, we remove the unphysical 
single unresolved limits of these four-particle tree-level antennae using appropriate products 
of two tree-level three-particle antennae, as was done in the massless case in [32,38]. 

In eq.(2.3), the four-parton antenna functions X^jj^-i depend on the original momenta 
Pi, kj, kk and ki while the reduced m-parton matrix element depend only on the redefined 
on-shell final-state momenta p^, . . . ,pl, . . .Pm+4 and on the rescaled initial-state momenta 
XiPi and p2- Thus, in order to obtain the integrated form of our subtraction terms, we 
factorise the full (m -|- 2)-particle phase space (with parton i in the initial state) in the 
following way [2]: 

d^m+2{p3, ■ ■ ■ ,Pm+4.;Pi,P2) = d$„i(?'3, • • • ,PL, ■ ■ ■ ,Pm+4.; XiPi, P2) 

X — —d^3{Pj,Pk,Pi;Pi,P2) , (2.4) 

with m?jj.i = mj + m| -|- rrii and 

2pj • q 

Using the phase space factorisation in eq.(2.4), the integrated form of the parts of the 
subtraction terms that involve four-parton antennae are given by 

\Mm+2?j'^'^ d^m I ^-^p^d^s{pj,Pk,Pi;Pi,q) (2-6) 

The integrated form of the four-parton initial- final antennae, denoted by X^jj^i, are 
thus obtained by integrating the four-parton antennae X^^j^^ over the 2^3 particle phase 
space d^3{pj,pk,pr,Pi,q) (with parton i in the initial state), analytically in d = 4 — 2e 
dimensions. More precisely, the integrated four-parton initial-final antenna is defined as: 



-^ijkli^i) = 7^-772 J d^3ipj,Pk,Pi;Pi,q)Xfjkl , (2.7) 



1 + m 

[C{e)] 

where 



C{e) = {AttY^. (2.8) 



-67 

8^ 

The kinematics of the reduced matrix clement appearing in eq.(2.3) depends on the 
masses of the mapped momenta and on the momentum fraction carried by the remapped 
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initial state parton, Xi. As will be shown in section 4, the integrated antennae also depend 
on Xi, in addition to their natural dependence on the mass of the hard radiator 

3. Antenna functions for double real radiation 

In this section we define the three initial-final four-parton massive antenna functions whose 
integrated form will be given in section 4. These antennae have been derived for the first 
time in [2,40], and their explicit expressions will be given in appendix A for completeness. 
They are present in the subtraction terms needed for the computation of the partonic 
processes qq — ?• QQq'q' and gg — > QQqq, which are part of the double real radiation 
contributions to heavy quark pair production in hadronic collisions. 

In general, an antenna function is characterised by its parton content, and its two 
radiators. These are the hard partons onto which the antenna collapses to in the unresolved 
limits. Accordingly, antenna functions are grouped into quark-antiquark, quark-gluon and 
gluon-gluon antennae and they are all derived from physical matrix elements related, in 
the final-final case, to the decay of a colourless particle into partons [43,44]. 

As stated before, while at NLO, only trcc-lcvcl threc-parton antennae are needed to 
capture the single unresolved behaviour of the real radiation matrix elements, at NNLO, 
four-parton tree-level antennae and one-loop three-parton antennae are also required. The 
former are needed to capture the double unresolved behaviour of double real matrix el- 
ements while the latter are used to capture the single unresoved behaviour of one-loop 
matrix elements. Restricting ourselves to the integration of the subtractions terms needed 
to evaluate double real contributions, we shall not discuss the one-loop antennae in the 
context of this paper. 

In contrast with the massless antenna functions, massive antennae have explicit mass 
terms and can be of two different natures: flavour-conserving and flavour-violating. The 
latter are derived from partonic processes with a flavour-violating vertex connecting ra- 
diators of different flavours, one of them being massless and the other massive. Massive 
three-parton flavour violating antennae have been derived in [39,40], and one specific quark- 
antiquark four-parton antenna of this kind has been derived in [2] . Since its integrated form 
will be derived in section 4, we shall here recall how it is defined. 

When considering the partonic process qq — >■ QQq'q' presented in [39] , to account for 
the singularities that arise when the massless final-state quark-antiquark pair is unresolved 
between a massive final state fermion and a massless incoming one, our subtraction terms 
need a massive initial-final fiavour-violating quark-antiquark B-type antenna denoted by 
Z?4(1q, 4q, 3g, 2q'). In this antenna, the hard radiators are a massive final state quark (1q), 
and a massless initial state quark (2^/) ^. 

This antenna is derived from the matrix element of the process V*q' — >■ Qqq, with a 
flavour violating vertex joining the colourless virtual particle V*, an initial state massless 
quark q' and a massive final state quark Q. The explicit expression of this antenna together 
with its infrared limits has been given in [2]. It is given in appendix A for completeness. 

^In the explicit expressions of the antenna functions, we shall label initial state partons with a hat, 
massless fermions as q and massive fermions as Q. 
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The other four-parton antennae whose integrated form is derived in section 4 are 
two kinds of quark-gluon antennae. Those are employed in our subtraction terms for the 
process gg — )■ QQqq [40], in order to reproduce the double unresolved limits of the double 
real radiation matrix element in which the final state qq pair becomes unresolved between 
the massive (anti) quark and an initial state gluon. These two quark-gluon antennae are 
obtained in their final-final forms from the process x — ggqq, with the massive gluino 
g playing the role of the massive (anti) quark of mass tuq. The initial- final antennae 
required here are then obtained by crossing the gluon to the initial state in these final-final 
ones. The full amplitude for the process x — ^ 99QQ contains leading and subleading colour 
pieces [43]. By squaring the leading colour piece, in which the qq pair is emitted between 
the gluino and the gluon in the colour chain, one-kind of quark-gluon antenna, namely the 
antenna, is obtained, while, squaring the subleading colour piece, in which the gluon is 
emitted between the qq pair, the antenna is obtained. The antenna £^4(1^,3^,4^,2^) 
accounts for the infrared limits associated to the emission of an unresolved qq pair between 
a massive (anti) quark and an initial state gluon, while £4(1^, 3g, 4^, 2g) is used to account 
for the triple collinear limits that involve a massless qq pair and an initial state gluon in 
those sub-leading colour amplitudes in which the gluon is placed between the quark and 
the antiquark in the colour chain. 

The explicit expressions of these two quark-gluon initial-final antenna functions de- 
noted by E2{lQ,3g,4:g,2g) aud E2ilQ,3g,ig,2g) havB been derived together with their 
infrared limits in [40] and are recalled in appendix A for completeness. 



4. Integration of initial-final antenna functions at NNLO 

In this section we describe our calculation of the integrated initial-final four-parton anten- 
nae denoted by B^, q^^, £g Qgg and £gQqq whose unintegrated forms, defined in section 3 
are denoted by 54(1q, 4g, 3g, 2g/), ^^(Iq, 3,, 4g, 2^), and £4(1^, 3g, 4g, 2^) respectively. 

4.1 General structure of the integrated antennae 

The initial-final antennae considered in this paper have a DIS-like 2 — >■ 3 kinematics 

9 + +P3 +P4 (4.1) 

with p\ = rriQ, P2 = P3 = p1 = and q^ < 0. It turns out to be convenient to parametrise 
this kinematics in terms of the following variables 



2p2 ■ q E^^ 



where E^^ = {pi + P'i + Pa)"^ = {q+P2)^- The more familiar variables Ecm and mq can be 
related to Q^, y and z through 

i7i2 y /n2 2 y^ ^2 /A o\ 

Ecm = Q > "mq = Q ■ (4-3) 

1 — y — z ^ 1 — y — z 
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The fact that the Ecm > itt-q imphes that < z < 1 and y > and, rccalhng that > 0, 
it follows from eq.(4.3) that y < 1 — z. Thus, in terms of our variables y and z the physical 
region is given by 

0<z<l, 0<y<l-z. (4.4) 

Following eq.(2.6) we integrate the initial-final four-parton antennae over the three 
particle phase space and express our results in the variables Q, y and z: Xf-f^i{Q'^,y,z). 
The integration is carried out following the standard technique of reduction to master 
integrals using integration-by-parts identities (IBP) [45,46]. 

We start by expressing our phase space integrals as cuts of two-loop four-point functions 
with two off-shell legs in forward scattering kinematics [47] , and then reduce these two-loop 
integrals using the Laporta algorithm [48] as implemented in FIRE [49]. 

In order to write the phase-space integrals of our four-parton antenna functions as 
two-loop integrals with cuts, we consider the following propagators: 

Di = pI- rriQ 

D2=pI 

D3=pI 

Di = (p3 + P4 - P2f = {q- Pl)^ 

D5 = {P3 +P4f 
Dq = {P3 -P2f 
Dj = {pA -P2f 

Ds = {pi +P3f - rriQ 

Dg = {P3 +P4- qf -mQ = {Pi - P2f - rriQ, (4.5) 
where Di, D2 and D3 are cut in the phase space integration. In the reduction procedure 

Q 0000 

we impose momentum conservation q+P2 = P1+P3+P4, we set pf = rriQ, P2 = P3 = Pt = ^ 
and q'^ = —Q^, and we discard all integrals in which the cut propagators Di, D2 or D3 are 
not in the denominator. 

After carrying out the reduction, we find that the NNLO integrated four-parton an- 
tennae B'^, Q^^, £g and £^ can be expressed in terms of four master integrals 

-f[o] = J d^3iPi,P2,P3;P2,q) (4.6) 
/[_8] = J d<^3iPi,P2,P3;P2,q)iiPi +P3f -mq) (4.7) 

-f[4] = j d$3(pi,y2,P3;y2,g) -p^)2 ('^■^) 
%,9] = / d$3(Pl,i?2,P3;P2,g)7 T2 2T- (4-9) 

All these master integrals contain multiplicative factors of the form y^^^ which regulate 
soft endpoint singularities in initial state convolution integrals. These factors ought to be 
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-^[0], -^[-8] I [4] I [4,9] 

Figure 1: The topologies and mass distributions of tlie four master integrals encountered in the 
calculation of the integrated NNLO four-parton antennae B^, q^^, £^ Qqq and £^ Qqq- Bold (thin) 
lines refer to massive (massless) scalar propagators. The double line in the external states represents 
the off-shell momentum g, with — ~Q^. The cut propagators are the ones intersected by the 
dashed line. 



left as such in the master integrals themselves. All other terms in the master integrals can 
be expanded in e. Explicitly the master integrals which we denote collectively by /a(y, z, e) 
take the form: 

/„(y,z,e) =?/'"-2^i?„(y,z,e). (4.10) 

The integer m is specific to each master integral and the function Ra{y, z, e) is regular as 
y — 7- and can be calculated as Laurent series in e. 

The integrated antennae collectively denoted by X{y,z,€) are linear combinations of 
these master integrals with coefficients containing poles in e as well as in y. After the 
masters have been inserted into the integrated antennae, these take the form 

X{y,z,e) = y-^-^'np,{y,z,e) (4.11) 

where TZx{y, z,e) is a regular function as y ^ 0. The e expansion of the singular factor 
is done in the form of distributions: 



\m 

-l—ne 



, __*M+J^(z!^f !,.„(,) (4.12) 

ne ^-^ mi 

m=0 

with 

It is worth noting that in the functions Raiy,z,e) and TZx{y^ z,e), which are regular 
as y — 7- 0, the limit z — > 0, corresponding to the massless limit, cannot be safely taken. 
In those functions, terms proportional to log(z) = log{mQ / E'^^) are present. These terms 
are expected to be there and correspond to the quasi-collinear kinematical configuration 
[17,39,42]. 

To evaluate the integrated antennae, we distinguish two regions depending on the value 
of y: a hard region where y 7^ and a soft region where y = 0. The highest order in e 
needed in the expansion of each master integral is determined by the e and y-dependent 
coefficient that multiplies the integral in the integrated antenna. In the soft region, since 
the expansion in distribution generates an additional 1 /e factor, the function i?„ is required 
one order higher in e than in the hard region. 
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4.2 Master integrals 

In this part of the section we present the results concerning the four masters required for 
the evaluation of the integrated antennae defined above. While and /[_§] have already 
been calculated in [42], /[4j and /[4 9] are presented here for the first time. 

4.2.1 The master integrals /[q] and /[_§] 

The master integrals /[g] -^[-8] do not contain any propagator involving the incoming 
momenta q,p2- Therefore, they are effectively 1^3 phase space integrals with incoming 
momentum q+p2 and a massive particle in the final state. These integrals were computed 
in [42] in the context of the evaluation of final-final three-parton antenna functions with 
massive final state fermions. They are given by the following all order expressions 



X2-Fi(l-e,2-2e,4-4e;l-z) (4.14) 



X2Fi(l-e,2-2e,5-4e;l-z), (4.15) 

with 

The hyper geometric functions in eqs.(4.14) and (4.15) can be expanded with the Mathema 
tica package HypExp [50], yielding ordinary harmonic polylogarithms (HPLs) [51] in the 
variable z. In the computation of the integrated antennae of this paper, we find that these 
two integrals are needed in the hard regions up to O(e^) while in the soft region, they are 
needed up to 0{e^). At this order (0(e^)), the harmonic polylogarithms (HPLs) in the 
variable z are needed up to weight 3. 

4.2.2 The master integrals /[4] and /[4 9] 

To evaluate the two remaining master integrals we use the method of differential equa- 
tions [52-58]. In order to obtain those we proceed as follows: We start by noting that 

^ 0^2 9 1-y-z-yzd 1 - 2y - z d 



dq^ dQ'^ 1 — z dx y dz 

d ,^ , d z(l — z) d 
P2fMj, — = (1 -y)^ - 



'dp2,^ dy y dz 

drriQ 1 — zdy'^^dz^ (^•-'-^) 



and invert this system to get 



2 5 19,1 9 , Q'^yz 



Q -^7^ = TiQt^a ^ ^P2l^a ^ 



dQ'^ 2 dqn 2 5p2/i 1 - y - ^ dm'^ 
d 1 d Q^z{l-z) d 



dy 1-y - z^^^ dp2^ ' (1 - y - z)^ dml 
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d 



dz {1 - z){l - y - z)^'^^ dp2^ ' (1 - y - z)"^ dml,' 



(4.18) 



Q 



We then apply the differential operators in eq.(4.18) to the integrals /[4] and /[4 gj and reduce 
the resulting integrals on the right hand side to master integrals. We thus obtain a set of 
first order partial differential equations for the master integrals which can be readily solved 
with standard techniques. The solution is found order by order in a Laurent expansion 
in e. This expansion involves harmonic poly logarithmic functions of one variable (HPL's), 
which can be functions of either y or z, or generalised harmonic polylogarithmic functions 
(GHPL's) [59,60] of variable z and arguments which can are taken from the list (0, 1, 1 — y). 
Products of these two type of functions are also found. 

The master integral /[4] is finite. For the integrated antennae required in this paper, 
it is needed in the hard region (y ^ 0) up to ©(e^) while in the soft region it is needed 
up to Oie'Y At this order (0(e^)) it will involve polylogarithmic functions, whose overall 
weight can go up to 4, as the 0{^) involves polylogarithmic functions of weight 1. This 
master integral can, however, be calculated in both regions up to order e^. In the hard 
region (y ^ 0) it is given by 

= \C{^)f (Q' + m%) |y-2^(l - y - z)'' | - (1 - y)(l - z)iJ(l; y) - yzH{0; z) 
+e (1 - y){l -z)(^- 2H{1; y)G{l -y;z)- 2G{1 - y, 0; z) - 2H{1; y)H{l- z) 
-2i/(0, 1; y) + ZH{1, 1; y)^ + yz (^2H{1; y)H{0; z) - 5H{0; z) + H{{), 0; z) 
-5H{1; y) - 2H{0, 1; z) + - 2{1 - y - z)H{l, 0; z) 

(1 - y){l -z)( - WH{1- y)G{l -y;z)- 4i/(0, 1; y)G{l - y- z) 



+6/f(l, 1; y)G{l -y;z)+ AH{1- y)G{l, l-y;z) + 4H{1; y)G{l - y, 0; z) 
+4H{1; y)G{l - y, 1; z) - 4H{1; y)G{l -y,l-y;z) + ^7r''G{l - y; z) 
-10G(1 - y, 0; z) + 4G(1, 1 - y, 0; z) + 2G(1 - y, 0, 0; z) + 4G(1 - y, 0, 1; z) 
-4G(1 - y, 1 - y, 0; z) - 10iJ(l; y)iJ(l; z) - 4H{0, 1; y)iJ(l; z) 
+6H{1, 1; y)H{l; z) - 4H{1- y)H{l, 1; z) - 19H{1; y) - 10/7(0, 1; y) 

+15i7(l, 1; y) - 4i7(0, 0, 1; y) + 6i7(0, 1, 1; y) + 2H{1, 0, 1; y) - 7/7(1, 1, 1; y) 

+yz [lOH{l- y)H{0; z) - 4H{1, 1; y)H{0; z) - 2H{1; y)H{0, 0; z) 

+4Hil; y)HiO, 1; z) - ^7r2//(l; y) + i (vr^ - 114) i/(0; z) 

+5H{0, 0; z) - 10/7(0, 1; z) - H{0, 0, 0; z) + 2H{0, 0, 1; z) - 4H{0, 1, 0; z) 

-4H{0, 1, 1; z) - 4H{1, 1, 0;z) + ^ (I2C3 + Stt^) j 
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+(1 - y - z) (^4H{1; y)H{l, 0; z) + \'^''H{1- z) - 10iJ(l, 0; z) 
+2H{1, 0, 0; z) - AH(l, 0, 1; z) - ^-'k''zH{1- y) 
+e3 (1 _ _ ^) _ (38 + H{l-y)G{l - y; z) 

-20i?(0, 1; y)G{l -y;z)+ 30H{1, 1; y)G{l - y; z) - 8H{0, 0, 1; y)G{l - y; z) 
+12H{Q, 1, 1; y)G{l -y;z)+ AH{1, 0, 1; y)G{l -y;z)- UH{1, 1, 1; y)G{l - y; z) 

+^ (Stt^ + I2C3) G{l-y-z) 

-^7r2G(l, \-y-z) + ^ (-114 + TT^) G(l - y, 0; z) 

+^^2^(1 _ y, 1 _ ^) + 20G(1, 1 - y, 0; z) + 10G(1 - 0, 0; z) 

+20G(1 - y, 0, 1; z) - 20G(1 - y, 1 - y, 0; z) - 8G(1, 1, 1 - y, 0; z) 

-4G(1, 1 - y, 0, 0; z) - 8G(1, 1 - y, 0, 1; z) + 8G(1, 1 - y, 1 - y, 0; z) 

-2G(1 - y, 0, 0, 0; z) - 4G(1 - y, 0, 0, 1; z) + 8G(1 - y, 0, 1, 0; z) 

-8G(1 - y, 0, 1, 1; z) - 8G(1 - y, 1, 1,0; z) + 8G(1 - y, 1, 1 - y, 0; z) 

+4G(1 - y, 1 - y, 0, 0; z) + 8G(1 - y, 1 - y, 0, 1; z) - 8G(1 - y, 1 - y, 1 - y, 0; z) 

+20G(1, 1 - y; z)H{l- y) + 20G(1 - y, 0; z)/i'(l; y) + 20G(1 - y, 1; z)i/(l; y) 

-20G(1 - y, 1 - y; z)/f(l; y) - 8G(1, 1, 1 - y; z)H{l- y) 

-8G(1, 1 - y, 0; z)H{l- y) - 8G(1, 1 - y, 1; z)H{l- y) + 8G(1, 1 - y, 1 - y; z)H{l- y) 

-4G(1 - y, 0, 0; z)if(l; y) - 8G(1 - y, 0, 1; z)H{l- y) - 8G(1 - y, 1, 1; z)H{l- y) 

+8G(1 - y, 1, 1 - y; z)i/(l; y) + 8G(1 - y, 1 - y, 0; z)F(l; y) 

+8G(1 - y, 1 - y, 1; z)i?(l; y) - 8G(1 - y, 1 - y, 1 - y; z)H{l- y) - 65H{1; y) 

-38H{1; y)H{l- z) + 8G(1, 1 - y; z)i7(0, 1; y) + 8G(1 - y, 1; z)i/(0, 1; y) 

-8G(1 - y, 1 - y; z)if(0, 1; y) - 20if(l; z)if(0, 1; y) 

(114 - 57r2) i?(0, 1; y) - 12G(1, 1 - y; z)Ji-(l, 1; y) 
-8G(1 - y, 0; z)H{\, 1; y) - 12G(1 - y, 1; ^)i7(l, 1; y) + 12G(1 - y, 1 - y; z)H{l, 1; y) 
+30iJ(l; z)iJ(l, 1; y) + 57i7(l, 1; y) - 20F(1; y)F(l, 1; z) - 8H{0, 1; y)i7(l, 1; z) 
+12H{1, 1; y)/f(l, 1; ^) - 8F(1; ^)F(0, 0, 1; y) - 20F(0, 0, 1; y) 
+12iJ(l; ^)F(0, 1, 1; y) + 30i7(0, 1, 1; y) + 4i7(l; z)H{l, 0, 1; y) 
+10F(1, 0, 1; y) - 14i7(l; z)H{l, 1, 1; y) - 35i?(l, 1, 1; y) - 8F(1; y)H{l, 1, 1; z) 
-8H{0, 0, 0, 1; y) + 12i7(0, 0, 1, 1; y) + 4H{0, 1, 0, 1; y) - 14F(0, 1, 1, 1; y) 

+4/f(l, 0, 0, 1; y) - 6iJ(l, 0, 1, 1; y) - 2H{1, 1, 0, 1; y) + 15iJ(l, 1, 1, 1; y) 
+yz (h{0, 0, 0, 0; z) + 4il(0, 1, 1, 0; z) + 4il(l, 1, 0, 0; 
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_1 (_ii4 + TT^) H{0; z)H{l; y) - \''^''H{1- y) 
+^'K''H{l;y)H{l-z) - 10i?(l; y^O, 0; z) 
+^ (114 - vr^) H{{), 0; z) + 201^(1; y)H{Q, 1; z) 



38 



137r^ 



F(0,l;z) -20/7(0; z)i/(l,l;y) 



+4i7(0, 0; z)H{l, 1; y) - 8/7(0, 1; z)H{l, 1; y) + ^t^''H{1, 1; 

+2i7(l; y)i7(0, 0, 0; z) - 5H{0, 0, 0; z) - 4/7(1; y)H{0, 0, 1; z) 

+1077(0, 0, 1; z) + 877(1; y)77(0, 1, 0; z) - 2077(0, 1, 0; z) + 877(1; y)77(0, 1, 1; z) 

-2077(0, 1, 1; z) + 877(1; y)77(l, 1, 0; z) - 2077(1, 1, 0; z) + 877(0; z)77(l, 1, 1; y) 

-277(0, 0, 0, 1; z) + 477(0, 0, 1, 0; z) + 477(0, 0, 1, 1; z) + 477(0, 1, 0, 0; z) 

-877(0, 1, 0, 1; z) - 877(0, 1, 1, 1; z) - 877(1, 1, 0, hz) + ^ (IMtt^ - tt"^ + SGOCa) 
+^77(0; z) (-390 + Stt^ + I6C3) 
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+{l-y-z)i^- 877(1, 0, 1, 0; z) + 877(1, 1, 1, 0; z) - n^H{l- j/)77(l; z) 
+2077(1; y)77(l, 0; z) - ^ (ll4 - vr^) 77(1, 0; z) 

-877(1, 0; z)77(l, 1; y) - 477(1; y)77(l, 0, 0; z) + 1077(1, 0, 0; z) + 877(1; y)77(l, 0, 1; z) 
-2077(1, 0, 1; z) - 277(1, 0, 0, 0; z) + 477(1, 0, 0, 1; z) - 877(1, 0, 1, 1; z) 

+ ^77(1; z) (Svr^ + I2C3) - ^vr277(l, 1; y) 



+ 



fvr2 + fC3)i^(l;t/) 



+0(0 ■ 



(4.19) 



The expression of 7[4] in the soft region is simply obtained by setting y = in eq.(4.19). 

The master integral 7[4 9] appears in the reduction of the integrated antennae into 
masters. It is finite, and a priori only its soft limit is needed up to order 0(e°). Up to this 
order, its expression in the hard region reads 



Im = [C{e)f[Q' + miy''^-y~'^{l-y) 



77(l;y)77(0;z) 



-77(1; y)G(l -y;z)- G(l - y, 0; z) + 77(0, 1; y) + 77(1, 1; y) - 77(1, 0; z) + 0(e) 

(4.20) 

With the soft limit of this expression being zero, it turns out that this master integral does 
not contribute in the result obtained for the integrated antennae derived in this paper. 
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4.3 Integrated antennae 

The three integrated antennae q^^, and f ^ q^- have a common multiphcative 

factor given by (Q^ + tjIq)^'^''. The full expression of these integrated antennae are too 
lenghty to be presented here. A separate Mathematica file containing them is attached 
with the arXiv submission of this paper. In this section we shall present the complete pole 
part of these integrated initial-final four-parton antennae. 

The leading pole parts of B^, and Sg Qqq are obtained in the kinematical configura- 
tion where the final-state qq pair becomes soft and simultaneously collinear to an initial- 
state parton. Those leading poles parts arc proportional to l/e'^5{y). 

Having no soft qq singularities, the integrated antenna £g Qqq is regular in the limit 
y — )■ and it does not contain any distributions. Its leading pole piece is proportional to 

To make our expressions more concise we define the following variable 



X = J,^^±- (4.21) 



which is equal to mq/Q. 

Using this notation, the pole parts of the integrated antennae are: 



1 



1 

+- 

e 



■, , , / 373 Stt^ I9z 25^2 \ ^ , , /19 z 
^^y^ (,-432 + ^-^ + 144;+^°^^^ (,36 + 6-12 

1 ^ , , 11 5y z ^, . f 1 z z'^\ 

-s"'*"' - 18 + 3I + T3I*fe) (6 + 8-24) ^> 

/IQ Z z'^\ 1 1 

^l^oivmUz) + (1 - I) HiO;y) - - ^ + |) m;y) - (i - |) i^(l;.) 
+O(e0)j, (4.22) 



^ ^-Ay) 



-2e 



X< - 



e3 12 



+- 



1 



^^^)(-72-T2 + 24 +6^°(^) 



1 



6 



5{y)Hil;z) 



+- 



1 



6 8 



2V 
24 



+ 



H{l;y) 
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+- 



4 



373 57r2 
432 ^ 



19z 25z' 



72 ^ 144 



+ Do{y) 



19 z 

36 ^ 6 ^ 12 



1 _ A 7y _ 29y2 



Xy'^ 4y^ z z^ 
^ ~3~ ~ 6 ~ 12 



y 



l-x 



^(l-j/)(l-^)V3^ 3 
^(^^G-4(l^ + ^ 



(l-j/)(l-z)2V3 

^) H{^- z) - S{y) ^^ + ^-^^H{l;z)- ^<5(y)iJ (1, 0; z) - ^<5(y)i7 (1, 1; z) 



+-Doiy)Hil;z) + { 1 



1/(1 - 
2(1 -y-z) 



1 



^ + ^)i?(i;y) 



1 _ A _ 2 
+ 6 ~ 2 ~ 3(1 -y) 3y 

5A_ 3y V 

"^1 3~ 2(1 -y-z) 



+ 



2 1 1 

3 ~ 2 + 6(l-y) 



y(i - y) 



+ 



l-z 



3y llA 



4 



2 - 13A 
6(1 -y) 
(2-3Ay 



1 



(1 - zY^ 



_ 7 _ 3y 7A 



+ 2Xy-y^ + Xy^+^^'-^^^ 



6(1 -y) 



2(l-y)(l-2^3 



H{0;z)+ -1 + 



5y 2y2 



3(1 -y) + 4 



1^(1; z) 



i/(l;y)ff(l;z) 



Hil;y)Gil-y;z)+{l 



Hil,0;y) 



+ 

+0(e° 



f) 



i^(i,i;y)- 1 



H{l,0;z)-[l-^]G{l-y,0;z) 



(4.23) 



2 \-2£ 



+ 



+ 



iiy y' 



3(1 -y) 12 3 



i-^ i?(i;y) 



1 

+- 

e 



20 ^ 16y ^ . 



13y2 1 



1 



+ 



y^ 

(i-y)(i-^) 

p lly 2y^ 2 
^lv3 6^3 3(1 -y) 

..2 



y-y 

l-y-z 

y 



H{l;y) + 



y 



+ 



(i-y)(i-^) 

..3 /4 



y- 



(l-y)(l-z)3V3 



1-y 

3y + 2Ay) + 
x))h{0;z) + 



H{0;y) + 



2 + A - 3y + 



1 , 7y 2y2 
3"^" 12 + ~ 



2y 



3(1 - y) 



y(i - y) 



y^ 



(l-y)(l-z 

y (^15 



1 - y - ^ 

(_l + 2A-2y + 2Ay) 



15y ^ 2y^ , 

(i-y)V6 + 
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+(2 - y)H{l; y)H{l; y) - (2 - y)H{l; y)G{l -y;z)+\^3- -f j H{1, 1; y) 
+(2 - y)H{l, 0; y) - (2 - y)H{l, 0; z) - {2 - y)G(l - y, 0; z) 
+O(e0)|. (4.24) 

5. Conclusions 

We have derived the integrated forms of the massive initial-final tree-level four-parton 
antennae appearing in the subtraction terms required to compute the partonic processes 
qq ttq'q' and gg — >■ ttqq contributing to the double real corrections to tt hadronic 
production. We found that those integrated antennae can be written as combinations 
of three master integrals from which one was unknown. This master integral is derived 
using differential equation methods and presented for the first time in this paper. These 
integrated antennae, related to the terms proportional to the colour factor Np in top 
quark pair production at NNLO represent the core result of this paper. With the integrated 
forms of these antenna functions now available, the integrated subtraction terms containing 
these antennae can be easily computed and combined with the double-virtual, one-parton 
integrated real-virtual and mass factorisation counterterms participating in the m-parton 
final state channel of an m-jet observable evaluated at NNLO. These integrated subtraction 
terms will contribute in a non-trivial way to the cancellation of the explicit poles present 
in this m-parton channel. As such, the results presented in this paper are essential for 
the application of the extended antenna formalism to compute hadronic observables with 
massive fermions at the NNLO level. These results, furthermore enable the computation 
of the Np contributions to the hadronic ti cross section evaluated at NNLO. 



6. Acknowledgments 

We would like to thank Werner Bernreuthcr for comments on our manuscript. Oliver 
Dekkers would like to thank the institute for theoretical physics at ETH Zurich, where most 
of this research project has been carried out, for its kind hospitality. This research was 
supported by the Swiss National Science Foundation (SNF) under contract PP00P2-139192 
and in part by the European Commission through the 'LHCPhenoNet' Initial Training Net- 
work PITN-GA-2010-264564', which are hereby acknowledged. Furthermore, the work of 
Oliver Dekkers was supported by the Deutsche Forschungsgemeinschaft (DFG), SFB/TR9. 

A. Four-parton antennae 

In this section we list the expressions of the unintegrated initial-final four-parton massive 
antennae defined in section 3. 

The flavour-violating initial-final four-parton B-type antenna is given by: 

Bl{'^QAq,^q',\') = \{ (S12S13 + -812514 + -813523 + •514S24) 

[Q'^+ml] y ^34^134 
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2 

■5345234 



1 / 

(S12S23 + S12S24 - S13S23 " S14S24) 



(2si2Si3Si4 + S13S14S24 + S13S14S23 - S13S24 - S14S23) 



„2 „2 
*34*134 

+ 22 (-2S12S23S24 + S13S23S24 + S14S23S24 - S13S24 - SMsia) 

^34*234 

H (2si2 + S12S23 + S12S24 - S12S13 - S12S14) 

S34S134S234 

H — 2 ( ~ ■513S24 - S14S23 + ■513S24 + S14S23 

■S34-5134S234 

-S13S14S23 - S13S14S24 + .S13S23S24 + S14S23S24 



-2si2Sl3S24 - 2S12S14S23) 



2^12 
S134S234 



+2ml + ^— + £11^) I + 0{e), (A.1) 

VSl34 S134S234 S34'Sl34 S34SI34 J) 

where = -{pi - P2 + Pz + Pa)'^, •§134 = si3 + si4 + S34, and S234 = -S23 - •§24 + S34- In 
this case the normaUsation is 

\M'2{7*q ^ Q)|' = 4(1 - e)(Q2 + ^^). (A.2) 



The initial-final quark-gluon four-parton E-type antenna is given by: 



{q 



2 + 

Si4 



S13S14 



+ 



Si4 



[-S13S23 + sis + su] 



Sl2S34'Sl34 
Sl4 



+ 



■5l2.S24.'5234 
4S11-SJ2''523 



S12S234 S12S34S24 
[2S14S23 + 2S14S24 + - sj^] 

[-S13S23 + sis + su] + [2si3 + 2si4 - S23] 



+ 



Su 



S34S134S234 S34S134 



S12S24 

+ 4si2Si4 + 4si4S23 + 4si4S24 



+8S12S24 + 45^2 + 4523*24 + 4s|4 



Si4 



S34S24S134 
Si4 



[2S12S14 + S14S23 " 2Si4 - S12S23 " S 



S34S134 



[-4S12S23 - 4S12S24 - 2Si2 " 4523*24 - 2S23 " 2534] 



Su _^ s 



14 



2 2 [-4S12S23 - 4si2S24 - 2Si2 - 4523*24 - 2^13 - 2s|4] 
S234 Si^sfs^ 



1 



Sl2'S34'5234 
1 



[S13S?4 + slsSU - S?3^23 + S?3 + «14«23 + sf^] 



S12S34 L 



S13SU - 3S13S23 - Sl3'S24 + 2Si3 + •514S23 " S14S24 
1 



+45^4 + S23S24 + sis 



S12S134 



S13S24 - S14S23 
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-2si4S24 - •S23S24 + S24 — [*13 + 2S14 - S23 - S24] 

■^^^ [2Sl4S23 + 8Su - 2si2S23 + 2s?2 + ^sla] 



S34Sl34'S234 
S23 



+ 



4si3Si4 - 8S12S13 + 4si3S23 + 4Si3 + 4si4S23 + 4si2S23 



+4s?2 



'23 



2 „2 



S-iaS 



34^234 



[-4si3Si4 + 4si2Si3 - 25^3 + 4si2Si4 - 25^4 - 25^2] 



^34 



4S13S12 - 4si3S23 + 4si3S24 - 2Si3 " 4si4Si2 - 4si4S23 



-4si4S24 + 4si2S23 " 4si2S24 - 2Si2 - 2S24 



+ - 



S34S134 



S12SU + 7S14S23 - S14S24 + 7Si4 + S12S23 + 7S12S24 

1 



+6Si2 + 2S23S24 + + 3S24 



+ 



S34S234 



6S12S13 - 2S13S23 + 4s 



13 



1 



+4si4S23 + 6Si4 + 5Si2 H [3si3 + 2si4 - 6si2 - 6S23 - •S24] 

-I S34 



+- 



.S23 



S24S234 
1 



[2si3Si4 - 2si2Si3 + 5^3 - 2si2Sl4 + S14 + •Sl2j 



S24SI34S234 



S12S13S14 - S13S14S23 - Sl3S?4 + S12S13S23 + 3Si2Sl4 
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[S12S13 - S12SU - Si4 + S 
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•S24S234 
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1 



2Sl2Sl3 - 2si3S23 - sj^ - 3si2Si4 - S14S23 + 4s 



14 



+ — [2si3 - 2S14 + S34] 
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'134 



[-2S12S23 - 2si2S24 - ^12 " 2523*24 - S23 " ■^24] 
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S134S234 

H [-3si3 + 6si4 + 3S23 - 2S24] + 2 

S134 

+mQm^ — h — — h 

S23 



+ 



S12 


^ 3 


Si4 






S24'5l34 


■S34 


S34S134 


2.S-23 


1 


4.S14 


•534S345 


S345 


S12S24 



S14S23 



+ 



[S24 - S13] 



+ 



S14S23 



S24 



Si4 



S12S24S34 S12S134 SI2S34SI34 



[S23 + 3S24] 



2s 



14 



Sl3 



S14S23 



2 [*23 + 524] H 

S12S34S134 S345S12 S24S345S12 



2si3S23 



+ 



2^23 



[Sl3 + S14] 



S34S345S12 S34S345S12 
2 r 

H 2~ ^13^23 + S14S23 " S13S24 + S14S24 

S12S|4 



-19- 



Sl2 



\ L2s12 - Sl4 + S23 + 3S24j H 



S24'5l34 S34 ■534'5l34 

[Sl2 + Sl4 + S23 + S24] H ^^^^ [Sl2 - S14] 



S34S345 



S134S345 

4S23 



[S23 + S14] + 



1 



S345 •534-S12 



S24S134S345 

[4si4 - 2S23 - 3S24] 



S34'Sl34S345 

2 

— [s\2 + 2si2S23 + 2S12S24 + ^23 + •S24 + 2523524] 

534^134 

H [2S23 - 524] + „ } _ [2si3 + 5S14S23 + 2S24 - S14S24] 



+ 



5125134 

2 



512S345134 



5 12 534 5 134 
2523 



[514523 + 514524] 



5125345345 



[-2S23 + ^13523 - 2S14523] 



2 [•^IS + •S14] 2~ [^13523 + 514523 " 513524 + 514524] 

512S34S345 512534 



2 2 

+ Z 2~ [*13 + 514 - 523 " 524] 



512 



'12 



2523 
5125345 



■-2 [513 - 2si3523 + 5i4 + 533 + 534 - 2S14S24] 

5i2534 



+0(e), 



S12 5125345134 



[S23 + 524] 



+ 



2ml ^ 
— [523 + 524] \ 

5125345134 J 



(A.3) 

where = -(pi+P3+P4-P2)^, rriy, = ^/Q^, 5134 = 513+514+534 and S234 = 534-523-524- 
The initial-final quark-gluon four-parton E -type antenna is given by: 



E4{lQ,3q,4q,2g) 



523524 



+2513534 + 2Si3 + 2514534 + 25^4 



- 2si25i3 - 2si25i4 - 2S12S34 + 25^2 



523 



5234*24 L 



2512513 - 2S12S14 



+5i2 + 2si3Si4 + Si3 + 5^4 



+2513514 + Si3 + Si4 



+ 



524 



5235234 



2512513 - 2si2Si4 + 5' 



12 



5234524 



-2513514 - 513523 - 2Si3 " ^14523 



4512513 + 2si2Si4 + 512523 " 25i2 
1 



5235234 



2si2Si3 + 4si25i4 



+512524 - 25i2 " 2513S14 - 513524 - 514S24 - 2s^4 



2524 ^ 4^234 ^ 2S23 ^ 4 ^ 4 



5235234 523524 5234524 523 524 



+ 0(e), 



(A.4) 

with Q^, m,^, 5i34 and 5234 given as above. Both E-type antennae are normalised to the 
tree-level two-parton matrix element 



(A.5) 
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